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Abstract 

We present a systematic study of the Gambier system, which in the continuous case is 
given by two Riccati equations in cascade. We derive the condition for its integrabihty 
and show that the generic Gambier system contains one free function. We also derive 
the Schlesinger transformations for this system which aUows in principle the systematic 
construction of the integrable cases. The above procedure is carried over to a discrete 
setting. We show thus how the discrete Gambier system can be expressed as a system of 
two homographic mappings in cascade. The integrable cases are obtained through the 
singularity confinement discrete integrabihty criterion. Finally the discrete Schlesinger 
transformations are also derived giving a handle to the construction of the integrable 
Gambier mapping. 



1. Introduction 

Linearisable equations play a very particular role in the domain of nonlinear sys- 
tems. Their nonlinearity is not a genuine one in the sense that elementary transfor- 
mations can reduce them to a linear system. Calogero coined the term C-integrability 
(which applies to both ordinary and partial differential equations) in order to describe 
this situation and to distinguish it from the more complicated situations of systems in- 
tegrable through 1ST methods (S-integrability) . The simplest example of a linearisable 
system is the Riccati equation [1] : 

x' = ax^ + bx + c (1-1) 

which can be transformed through the Cole-Hopf transformation x = —u'j (au) to the 
linear second-order equation: 

u" = {b + a'/a)u' - cau (1.2) 

Higher degree, first-order linearisable equations also exist [1] . Already limiting ourselves 
to equations of binomial type we find 

{x'f = {x-af{x-ki){x-k2) (1.3) 

where a is a function of the independent variable and /ci,/c2 are constants. Putting 
= {x — ki)/{x — equation (1.3) reduces to the Riccati: 

u' = ±]-{ki - a - {k2 - a)u^) (1.4) 

When we consider second-order equations, the situation becomes immediately richer, 
in the sense that there exist several types of linearisable equation [1]. The first (and 
simplest) one is the equation obtained if one computes the derivative of both sides of 
the Riccati (1.1). It is usually given in canonical form and reads: 

x" = -2xx' + hx' + h'x (1.5) 

In the same spirit, one can compute the derivative of (1.1) after having divided both 
sides of the equation by x. We find thus the equation (again given in canonical form): 

x" = — Jr (ax --\x' + a!x^ c' (1.6) 

X \ X/ 

Another linearisable equation does exist which can be linearized by a Cole-Hopf trans- 
formation just like the Riccati. Its form is: 

x" = -Zxx' - x^ + q{x' + x^) (1.7) 

Indeed, putting x = u' ju we can reduce it to the linear third order equation: 



(1.8) 



But the most interesting linearisable equation discovered at second-order is the one 
obtained by Gambier in his classification of second-order equations having the Painleve 
property [2]. The Gambier equation is usually given as a second-order equation for a 
single variable: 

n — 1 x'"^ n + 2 , , n — 2 x' o . , , , o 

X = h a XX + ox a x + [a — ao)x 

n X n n X n (^ q\ 

2aa\ , ^^'^^ 



+ [cn )x — ba 

V n J nx 

where n is integer, a, b, c are functions of the independent variable and a is equal to 
or 1. Its structure becomes clearer when one writes it as a system. We have in this case 
a system of two Riccati's: 

y' = -y^ + by + c (1.10a) 

x' = ax^ + nxy + a (1.106) 

Although the Gambier equation is always linearisable this does not mean it is always 
integrable. Indeed, in a system such as (1.10) of two equations in cascade we can always 
solve the first equation for y, obtain y{t) and inject it into the second equation. The 
latter can always be written as a linear, second-order differential equation for x. So, 
in principle, the problem can always be solved formally. The difficulty comes when 
one wishes to actually compute x, in terms of contour integrals, while y has an analytic 
structure that interferes badly with that of x. This is where the Painleve property comes 
into play. If we require that the system possess the Painleve property the integration can 
be performed and we can indeed obtain the solution for x{t) over the complex t-plane. 
Thus, the integrability of the Gambier equation will be closely related to its singularity 
structure. 

The fact that the integrability of the Gambier equation is related to the Painleve 
property allows us to obtain another interesting result. Just as in the case of the Painleve 
equations, it is possible to introduce transformations relating the solution of an equation 
with parameter n to one with parameter n + 1 or n + 2 (two different transformations 
do exist). Thus the Gambier system possesses Schlesinger transformations [3]. 

Another feature of the Gambier equation is that it can be integrably discretized 
[4,8]. Our approach follows closely the spirit of Gambier based on coupled Riccati equa- 
tions. In perfect analogy to the continuous case it is possible to introduce a system of 
two homographic mappings in cascade which represent the discretization of the Gambier 
system. The general form of this system is the following: 

ayn + b n 1 1 \ 

Vn+i = —7 (1.11a) 

cyn + d 

a^n+i = 7 —T\ —J (1.116) 

(eyn + QXn + 'nyn + 

The form (1.11) can be simplified through homographic transformations and the in- 
tegrable cases can be obtained through the application of the singularity confinement 



discrete integr ability criterion. As in the continuous case, although (1.11) is always 
linearisable it is not automatically integrable. The difficulty arises when one tries to 
compute Xn in terms of matrix products the elements of which contain y^- Some of 
these matrices are singular in such a way that degrees of freedom are irretrievably lost 
when yn has the wrong properties. Singularity confinement precisely means that these 
degrees of freedom are in fact recovered at some later stage. 

Quite expectedly, the Gambler mapping has Schlesinger transformations. Their 
derivation is based on a study of the singularities of the mapping (1-11) and their 
confinement. Thus the parallel between the discrete and continuous cases is perfect. 

In what follows, we shall present the singularity analysis of the Gambler equation 
and derive its integrable cases. Next, we study the discrete case and use the singularity 
structure in order to derive the Gambler mapping. Finally, we present the Schlesinger 
transformations of both the continuous and discrete Gambler systems which allow us, 
starting from some elementary case (n=0 or n=l), to construct recursively the Gambler 
systems for higher n's. 

2. The continuous Gambier equation 

The Gambier equation is given as a system of two Riccati equations in cascade. 
This means that we start with a first Riccati for some variable y 

y' = -y^ + hy + c (2.1) 

and then couple its solution to a second Riccati by making the coefficients of the latter 
depend explicitly on y: 

x' = ax"^ + nxy + a. (2-2) 

The precise form of the coupling introduced in (2.2) is due to integrability requirements. 
In fact, the application of singularity analysis shows that the Gambier system cannot be 
integrable unless the coefficient of the xy term in (2.2) is an integer n. This is not the 
only integrability requirement. Depending on the value of n one can find constraints on 
the a, 6, c, a (where the latter is traditionnally taken to be constant 1 or 0) which are 
necessary for integrability. 

The common lore [1] is that out of the functions a, 6, c two are free. This turns 
out not to be the case. The reason for this is that the system (2.1-2) is not exactly 
canonical i.e. we have not used all possible transformations in order to reduce its form. 
We introduce a change of independent variable from t to T through dt = gdT where g is 
given by = ^23^' ^ S^^S^ through x = gX and also Y = gy — The net result 
is that system (2.1-2) reduces to one where 6=0 while a remains equal to or 1. It 
is clear from the equations above that n must be different from 2. On the other hand 
when n = 2 the integrability conditions, if cr = 1, is precisely b = 0. So we can always 
take b = 0. (As a matter of fact in the case u = an additional gauge freedom allows 
us to take both b and c to zero for all n, even for n = 2). Thus the Gambier system can 
be written in full generality 

y' = -y^ + c (2.3a) 



x' = ax^ + nxy + cr. 



(2.36) 



One further remark is in order here. The system (2.3) retains its form under the 
transformation x ^ Xjx. In this case n — > — n and a and —a are exchanged. Thus in 
some cases it will be interesting to consider a Gambler system where a is not constant 
but rather a function of t. Still, it is possible to show that we can always reduce this 
case to one where cr = 1, while preserving the form of (2.3a) i.e. 6 = 0. To this end 
wc introduce the change of variables dt = hdT, x = gX and Y = hy — with 
h = a^^^'^~'^\ g = o""'/'^"'"^-'. With these transformations system (2.3) reduces to one 
with cr = 1 and 6 = 0. (In the special case n = 2, with 6 = 0, integr ability implies 
cr =constant, whereupon its value can always be reduced to 1). 

In order to study the movable singularities of the coupled Riccati system we start 
from the observation that from (2.3a) the dominant behaviour of y can only be j/ ~ 
l/{t — to). The next terms in the expansion of y can be easily obtained, and involve 
the function c and its derivatives. In order to study the structure of the singularities 
of (2.3b), we first remark that since the latter is a Riccati, its movable singularities are 
poles. However, (2.3b) also has singularities that are due to the singular behaviour of 
the coefficients of the r.h.s of (2.3b), namely y. Now, the locations of the singularities of 
the coefficients are 'fixed' as far as (2.3b) is concerned. However, from the point of view 
of the full system(2.3), these singularities are movable and thus should be studied. The 
'fixed' character refiects itself in the fact -1 is not a resonance. (The terms 'resonance' 
is used here following the ARS terminology [6] and means the order, in the expansion, 
where a free coefficient enters. A resonance -1 is related to the arbitrariness of the 
location of the singularity, and is thus absent when the location of the singularity is 
determined from the 'outside' rather than by the initial conditions). Because of the 
pole in y, X has a singular expansion with a resonance different from -1 which may 
introduce a compatibility condition to be satisfied. 

We consider below the case of the full Riccati (2.3b) with a 7^ 0, cr = 1. (The 
analysis of the case acr = was given in [4]). As we explained above, only the singularity 
due to y can lead to trouble. Rewriting (2.3b) as x' /x — ax + ny + a/x for y = 
l/{t — to) + . . . we remark that unless n = ±1 a behaviour of the form a; ~ (t — to)" 
is impossible when aa ^ 0. For n = 1, a logarithmic leading behaviour will be present 
for cr 0. (Note that the condition cr = is sufficient for the absence of a critical 
singularity for n = 1 irrespective of the value of a. Similarly, in a dual way, for n = — 1 
the necessary and sufficient condition for the absence of a critical singularity is a = 0, 
whether cr is 1 or 0). 

Next we assume n 7^ ±1, in which case it suffices to study the singularities x 
X{t — to), (A = cr/(l — n)) and x ^ ji/{t — to), {fi = — (n + l)/a)- The first singular 
behaviour {x ^ X{t — to)) has a resonance at n — 1, which is negative for n < 1 and thus 
docs not introduce any further condition. For n > 1, the resonance condition can be 
studied at least for the first few values of n. (In fact, cr = suffices for the resonance 
condition to be satisfied even for a ^ 0). In the particular case n = 2, we have already 



mentioned that the integrabihty condition for cr = 1 is precisely 6 = 0. For higher n 
(and cr = 1) we find the further possibihties: 

n = 3 2c-a = (2.5) 

n = 4 3c' - a' = 

The second singular behaviour, x ji/{t — to), has a resonance at —1 — n. Thus 
for n > this resonance is negative and does not introduce any further condition, while 
for n < a compatibility condition must be satisfied. (We find that for every case 
n < 0, a = is a sufficient condition for the absence of critical singularity. This is not 
in the least astonishing given the duality of a and a). On the other hand if we demand 
a 7^ then a different resonance condition is obtained, at each value of n. For n = — 1, 
whenever a 7^ 0, a logarithmic singularity of the form (t — to)~^log{t — to) appears 
irrespective of the value of a. For n < —1, we find: 

n=-2 a' = 



n=-3 2ac-aV-2a" = (2.6) 

n = -4 2ac' + a' {4c - 2aa) - a'" = 

The integrabihty condition for higher values of n can be obtained through the use of 
computer algebra. 



3. The discrete Gambier equation 

The discretisation of the Gambier equation is based on the idea of two Riccati 
equations in cascade. The discrete form of the first is simply: 

- ay + b 

y = — —7 (3.1) 

cy + d 

where y = y^ and y = j/n+i- The second equation which contains the coupling can be 
discretised in several, not necessarily equivalent, ways. In [5] we have considered the 
generic coupling of the form: 



axxy + (3xx + ^xy -\- 5x -\- exy + (x + rjy + 6 = (3.2) 

Implementing a homographic transformation on x and y we can generically bring (3.2) 
under the form: 

XX + jxy — exy — 9 = (3.3) 
. A choice of different transformations can bring (3.2) also to the form 



X — X = —fxx + {gx + hx)y + k 



(3.4) 



Note that (3.2) contains an 'additive' type coupling xx + 5x + Cx + r]y + 9 = for special 
values of its parameters, but the generic form (3.3) is that of a 'multiplicative' coupling 
where 7, e do not vanish. Solving (3.3) for x we obtained the second equation of the 
discrete Gambler system in the form: 

X = . (3.5) 

x + ^y 

Clearly, a scaling freedom remains in equation (3.5). We can use it in order to bring it 
to the final form: 

Eliminating y and y from (3.1), (3.6) and its upshift, we can obtain a 3-point mapping 
for X alone but the analysis is clearer if we deal with both y and x. 

In what follows we shall present a different derivation based on the study of sin- 
gularities of the system. As we explained in [3] we can use the homographic freedom in 
order to bring the mapping for y to the form: 

- y+c , . 

instead of (3.1) where c is a function of n. Next, we turn to the equation for x. This 
equation is homographic in x. However we require that when y takes the value 0, the 
resulting value of be 00. Thus the denominator must be proportional to y, and since 
we can freely translate x, we can reduce its form to just xy. The remaining overall gauge 
factor is chosen so as to put the coefficient of xy of the numerator to unity resulting to 
the following mapping: 

_^xiy-r) + qiy-s)^ (3.8) 
xy 

The system (3.7-8) is a discrete form of the Gambler system. In order to study the 
confinement of the singularity induced by y = we introduce the auxiliary quantity V-'n 
which is the A^'th iterate of y = in equation (3.7), N times downshifted. Thus •0o = 0, 

c+c 

V'l = c, ip2 = ^qij) etc... The confinement requirement is that after N steps x becomes 
in such a way as to lead to 0/0 at the next step. Thus the mapping (3.8) has in fact 
the form: 

_^ x{y-r) +qiy-ipN) gx 
xy 

Thus when at some step N we have y — ipN and a; = 0, on the view of (3.9) x will 
then be indeterminate of the form 0/0. However it turns out that in fact this value is 
well-determined and finite. Let us take a closer look at the conditions for confinement. 
The generic patterns for x and y are: 

y-{ Vi "02 • • • } 

X : { free 00 ^ ... free }. 



At A?" = 1 it is clearly impossible to confine with a form (3.9) since we do not have 
enough steps. In this case the only integrable form of the a;-equation is a linear one. 
This case must be studed separately. For an arbitrary N, the general form of the linear 
x-mapping can be obtained using confinement arguments in a way similar to what we 
did for the generic, nonlinear, case. We obtain 

y 

where g is free. 

The first genuinely confining case of the form (3.9) is A" = 2. From the requirement 
^ = we have r = V'l and q free: this is indeed the only integrability condition. For 
higher A's we can similarly obtain the confinement condition which takes the form of 
an equation for r in terms of q. 

At this point it is natural to ask whether the mapping (3.7)- (3.9) does indeed 
correspond to the Gambler equation (2.3). In order to do this we construct its continuous 
limit. We first introduce: 

c = e^D 

eD (3-11) 



y 



Y + H 

with H ^ D'/{2D) and obtain the continuous limit of (3.7) for e ^ 0. We find as 
expected 

Y' = -Y^ + C (3.12) 

(i.e. eq. (2.3a)) where C = D — ^ + Using (3.7) and (3.11) we can also compute 

V'AT and we find at lowest order: 

^^ = e2*jv with Wat ^ N{D - e ^^^ D') + e^^jv (3.13) 

where $jv is an explicit function of D depending on A. 
Next we turn to the equation for x and introduce: 

r = e^R 

1 e RD' 
"=2 + 2X-^4^ ^3.14) 

and for the continuous limit of the form (2.3b) to exist in canonical form (i.e. b=0, 
cr=l) we find that we must have 

ND , ,ND' 
i?«,___e(iV + 2)-^. (3.15) 

This leads to the equation for x: 



X' = AX^ + NXY + 1 



(3.16) 



with A = ^{N/A + V)%T — 4(5. Moreover the confinement constraint imphes 

a differential relation between D and Q which depends on N. We can verify explicitly 
in the first few cases that this is indeed the integrability constraint obtained in the 
continuous case. In the case N = 2, the condition for confinement is 

f = c. (3.17) 

Using the approximate expansion in e (3.15) which gives r up to the third order in e and 
putting N = 2, we get that equation (3.17) is automaticaly satisfied for the two first 
nonzero orders in e (e^ and e^). This reflects the fact that for the continuous Gambler 
equation, in the case N = 2, the only integrability condition is 6 = 0. Thus no condition 
has to be imposed on D and Q. 

The confinement condition for A'^ = 3 is obtained by imposing x = when y = 0. 
This condition reads: 

(c + c- (c+ 1)V^) = (f-c)(c + c-f(c+l)). (3.18) 

Putting N = 3, we find that equation (3.18) is identically satisfied for the first two 
nonzero orders in e (e^ and e^). To get the integrability condition, we must calculate 
(3.18) at order e^. To do this, we need the value of $3 in (3.13). We easily find that 
$3 = 7D" —8D^ . Satisfying (3.18) in gives a relation giving explicitly Q in terms of -D. 
Implementing this relation, we find that A — 2C which is the continuous integrability 
condition for A'^ = 3. 



4. SCHLESINGER TRANSFORMATIONS FOR THE CONTINUOUS GAMBIER EQUATION 

The theory of auto-Backlund transformations of Painleve equations is well estab- 
lished. As was shown in [7] the general form of auto-Backlund transformations for most 
Painleve equations is of the form: 

ax' + dx'^ + jx + d , , 

X = ■ ■ . (4.1) 

ex' + (^x'^ + rjx + 9 

In the case of the Gambler equation considered as a coupled system of two Riccati's it 
is more convenient to look for an auto-Backlund of the form: 

^ ^ axy + Px + ^y + 5 

{(y + rj)iex + k) ■ ^ ■ ^ 

with a factorized denominator, with hindsight from the discrete case. We require that 
the equation satisfied by x do not comprise terms nonlinear in y. We examine first the 
case C 7^ ^iid reach easily the conclusion that there exists no solution. So we take 
C = 0, 7/ = 1 which implies that a and 7 do not both vanish (otherwise (4.2) would have 
been independent of y). We find in this case a = and thus the general form of the 
auto-Backlund can be written as: 

_ px + ^y + 5 



From (4.3) we can obtain the two possible forms of the Gambier system auto-Backlund: 



X = (3x + ^y + 5 (4.4) 

x=^-^±jy±l. (4.5) 

As we shall see in what follows both forms lead to Schlesinger transformations. 

Let us first work with form (4.4). Our approach is straightforward. We assume 
(4.4) and require that x satisfy an equation of the form (2.3b) while y is always the 
same solution of (2.3a). The calculation is easily performed leading to: 

x = iy-\ — -xH , (4.6) 

n + 1 n 

where 7 satisfies: 

^ = ^-. (4.7) 

7 n + 2 o ^ ' 

Here we have assumed a ^ 0; otherwise x does not depend on x and (4.6) does not 
define a Schlesinger. The parameters of the equation satisfied by x are given (in obvious 
notations) by: 

n + n + 2 = (4.8a) 

a = ^ (4.86) 

7 

and 

CT = 7 cH \ zn^). (4.8c) 

'V n+1 n + 2 a {n+2f a? ) ^ ' 

Thus (4.6) is indeed a Schlesinger transformation since it takes us from a Gambier 
system with parameter n to one with parameter h = —n — 2. It suffices now to invert x 
in order to obtain an equation with parameter N = n + 2. Expressions (4.6) and (4.8) 
can be written in a more symmetric way: 

a' 

ax - ax = (n+ l)(y - —) (4.9) 

na 

and 

h+l = -(n+1) 
.a' a' 

^a=^a (4-10) 

/2- 



/ -, N / 1 / a X / 1 a ^\ 
aa - aa = {n+ l)[c - -{ — ) +^—;r] 
\ a' a^ / 



The inverse transformation can be easily obtained if we introduce 7 such that 
a7 = — (n + 1) = —07. We thus find 

x = y^+z:^x+%- (4.11) 
n + 1 n 



and the relations (4.10) are still valid. 



Iterating the Schlesinger transformations one can construct the integrable Gambier 
systems for higher n's and obtain by construction the functions which appear in them. 
However it may happen that when we implement the Schlesinger we find a = 0. If we 
invert x we get a system with N = —n = n + 2 but ^ = for which one cannot iterate 
the Schlesinger. 

Let us give example of the application of this Schlesinger transformation. Let us 
start from n = 0, in which case we find n = —2 and, after inversion, N = 2. For n = 
we start from a = — 1 and cr = or 1 (always possible through the appropriate changes 
of variable). This leads to a = —1, a = —c + a and the Schlesinger reads: x = —y + x. 
Next we invert x and have X — l/{x — y). We find thus that the Schlesinger takes us 
from 

y'r^M' (4.12) 

to the system 

y' = -y^ + c , . 

X' = + 2Xy + S ^^-^^^ 

with ^4 = — cr, E = l. In the particular case n = 2, a change of variables exists which 

allows us to put A = —\ (unless A = 0), without introducing h in the equation for y, 

while keeping S = 1 and changing only the value of c. Thus the generic case of the 

Gambier equation for n = 2 can be written with A = —\. Eliminating y between the 

two equations we find: 

a;" = 1^ - 2xx' - ^ - ^ + (2c + \)x. (4.14) 

This is the generic form of the n = 2, Gambier equation [1] and it contains just one free 
function. 

We turn now to the second Schlesinger transformation corresponding to the form 
(4.5). As we shall show, a Schlesinger transformation of this form does indeed exist and 
corresponds to changes in n with An = 1. Let us start from the basic equations (2.3). 
Next we ask that x defined by (4.5) indeed satisfy a system like (2.3). We find thus 
that and k = —xq and 7 must be given by: 

— = yo + ^— j- (4.15) 
7 n + 1 

where j/o is a solution of the Riccati (2.3a) and a solution xq of (2.3b), obtained with y 
replaced by j/o- We introduce the quantities xq = a = — In this case (4.15) 
becomes: 

i . 2axo 2x0X0 

— = yo + -^—T =yo + , (4.16) 

7 n + 1 7 

where 

n + n + 1 = 0. (4.17) 
We have thus, starting from a generic solution x, y of (2.3) for some n, the Schlesinger: 

x = xo + ^-^^^ (4.18) 

X — Xo 



where x is indeed a solution of (2.3) for h = 



—n — 1 for the same y 



x' = ax^ + nxy + a (4.19) 
where a has been defined as —uxq/j and 

c^= ^-r «' + «^^o^-r + ayo(n + 2) . (4.20) 
n+lV n+1 / 

Note that xq is a solution of the same equation with y replaced by j/o- As in the previous 
case if we invert x wc obtain an equation corresponding to = n + 1. 

It is worth pointing out here that the Schlesinger transformation corresponding to 
An = 2 was known to Gambler himself. As a matter of fact when faced with the problem 
of determining the functions appearing in his equation so as to satisfy the integrability 
requirement, Gambler proposed a recursive method which is essentially the Schlesinger 
An = 2. On the other hand the Schlesinger An = 1 is quite new. 

5. Schlesinger transformations for the discrete Gambier equation 

Once the singularity pattern of the Gambier mapping is established we can use it in 
order to construct the Schlesinger transformation. Let us first look for a transformation 
that corresponds to AA^ = 2. The idea is that given the A^-steps singularity pattern 
of the equation for x we introduce a variable w with + 2 singularity steps where we 
enter the singularity one step before x and exit it one step later. The general form of 
the Schlesinger transformation, which defines tu, is: 

w = X , (5.1) 

y 

where X is homographic in x. The presence of the y and y — i/jn+i terms is clear: they 
ensure that w becomes infinite one step before x, and vanishes one step after x. Next 
we turn to the determination of X. Since X is homographic in x we can rewrite (5.1) 
as: 

w = — . (5.2) 

y '^x + d 

Our requirement is that w becomes infinite when y — for every value of x. This 
statement must be qualified. The numerator ax + P will vanish for some x (namely 
X = —(3 /a) so this value of x must be the only one which should not occur in the 
confined singularity. Indeed there is a unique value of x where instead of being confined, 
the singularity extends to infinity in both directions of the independent variable n, while 
the only nonsingular values of the dependent variable occur in a finite range. The value 
of X such that x is finite and free even though y is zero is such that the numerator 
—xr — qi(^N ofx vanishes. For this value of x, the values of the dependent variable are 
fixed for n < and n > N + 1 and the value can be considered as 'forbidden'. Thus 
ax + /3 = xr + gViv up to a multiplicative constant. Similarly when y — i^N+i, w must 
vanish. Thus ^x + S must not be zero except for the unique value of x that does not 



occur in the confined singularity. Note that y = V'iv+i means y = i/^^ and the only- 
value of X that comes from a nonzero x in that case is a; = (V'^ — r)/i(^j^- In that case 
the values of the dependent variable are fixed for n > and n < —N — 1. This value of 
X being 'forbidden', + S must be proportional to V'^^; — — r). We now have the 
first form of the Schlesinger: 

xr + qipN y-i^N+i „x 

w = ■ ■ — (5.3) 

y + r-'ip^ 

where the proportionality constant has been taken equal to 1 (but any other value 
would have been equally acceptable). Here w effectively depends on x unless r(r — 
'ipj^) = Qipjyi^N- But in this case the mapping (3.8) is in fact linear in the variable 
^ = (x — 1 + r/ifjj^)~^. This case is the analog of the case a = in the continuous 
case where the Schlesinger does not exist. Let us give an application of the Schlesinger 
transformation by obtaining the N — 2 equation starting from N = 0. We have always 
the equation for y which reads: 

and V'o = 0, V'l = Q. For N = the equation for x reads: 

_ x{y-r) + qy x + q 

X = = (5.5) 

xy X 

since for integrability r = and indeed N = means that the x equation does not 
depend on y. We introduce the Schlesinger: 

y-c , . 

w = X (5-d) 

y 

Using (5.5) and (5.6) to eliminate x we obtain the equation for w: 

- .yw + q{y-c) 

w=(l-c) — -. r -. (5.7) 

This equation is of the form (3.8) but not quite canonical. We can transform it to 
canonical form simply by introducing y instead of y because indeed w is infinite one 
step before x, so w = oo means x = oo i.e. y — 0. We obtain thus: 

- ^ w{y - c) + q{l + c){y - ^2) .5 
yw 

with t(j2 = (c + c)/(l + c) which coupled to (5.4) is indeed a = 2 Gambler mapping. 

As we pointed out in section 3, necessitates a special treatment. The Schlesinger 
transformation is again given by: 



^ ^^ y-^N+i 

y 



(5.9) 



and arguments similar to those of the nonhnear case allow us to determine the form of 
the homographic object X leading to: 

W= ; . (5.10) 

y xipj^ - g 

Thus one can also perform a Schlesinger in the linear case. This is not in disagrement 
with the continuous case. It is, in fact, the analog of the case where (7 = but a 7^ 
(which is linear in 1/x) for which the Schlesinger can be performed. The analog of the 
case o" = and a = is the situation when g = kipN with constant k in which case the 
mapping rewrites ^ = ^{y — 'i(jN)/y with ^ = x — k. Then w does not depend on ^ (or 
x) and (4.20) does not define a Schlesinger in analogy to the case r(r — ip^) — qi^j^i^N 
in the nonlinear case. 

Finally, we examine the possibility of the existence of a AN = 1 Schlesinger. In 
this case, the structure of the transformation will be obtained by asking that the N +1 
case enter the singularity one step before the N case but exit at the same point. The 
general structure is thus: 

y x-i 

where r] and ^ must be determined. We do this by requiring that the equation for 
w contain no coefficients nonlinear in y. As a result we find that r] must satisfy the 
equation (3.7) for y: 

^ ,5.1.) 



and ^ the equation (3.9) for x with -q instead of y: 

- r) + g(?7 - V^iv) 



(5.13) 



We remark here the perfect parallel to the continuous case (and as we pointed out the 
discrete case led the investigation back to the continuous one). Let us point out here 
that the w obtained through (4.26) does not lead to w = at the exit of the singularity 
(i.e. when x = 0, j/ = '4>n) and a translation is needed. One has in principle to define a 
new variable 

q 

u) = w - w{x = 0,y = i/jn) = w + -{i/jn - v)- 

Finally we derive the AN = 1 Schlesinger for the case of a linear mapping (5.9). 
We start from: 

w = (5.14) 

y x-i 

and again require for w an equation with coefficients linear in y. We find that rj must 
again be a solution of the equation for y i.e. it must satisfy (5.12) and moreover ^ is a 
solution of (5.13) with y = rj: 

-^^av-M + g^ (5.15) 
V 



Thus the hst of the Schlesinger transformations of the Gambier mapping is complete. 
6. Conclusion 

In this work we have reviewed the Gambier system in both its continuous and 
discrete forms. We have shown that the singularity analysis can be used in order to 
obtain the integrable subcases of this equation and moreover we have derived their 
Schlesinger transformations. Several open questions appear at this point. The Gambier 
equation is the generic linearisable second-order equation of first degree. If one relaxes 
this last constraint, one can already obtain further linearisable equations. Cosgrove and 
Scoufis [8] have presented two such examples: 

{x"f = {ax' + a'x + c'fx' (6.1) 

{x")'^ = A{x){ci{tx' — x) + C2x'{tx' — x) + cs{x')'^ + C4{tx' — x) + c^x' + Cg). (6.2) 

It would be interesting to study closely the properties of these equations and in particular 
their discretization. 

Moving to higher orders one can wonder how the Gambier approach can be ex- 
tended. In [9] we have presented a first exploration of this problem, in both continuous 
and discrete settings, for third-order systems. 

Finally, the general problem of n-th-order linearizable systems is far from being 
solved. Only a special class of such systems is known, based on projective contructions. 
Recently, we have obtained the discretization of these projective Riccati systems [10]. 
Clearly, this is not the last word as far as linearisability is concerned. 
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